Double dot chain as a macroscopic quantum bit 
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We consider an array of N quantum dot pairs interacting via Coulomb interaction between ad- 
jacent dots and hopping inside each pair. We show that at the first order in the ratio of hopping 
and interaction amplitudes, the array maps in an effective two level system with energy separation 
becoming exponentially small in the macroscopic (large N) limit. Decoherence at zero temperature 
is studied in the limit of weak coupling with phonons. In this case the macroscopic limit is robust 
with respect to decoherence. Some possible applications in quantum information processing are 
discussed. 
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I. INTRODUCTION 

In recent years a lot of attention has been devoted to 
the existence of quantum superposition states in macro- 
scopic systems. The first suggestion to understand this 
phenomenon is due to Schrodinger Q who introduced the 
paradox of the cat in quantum superposition between live 
and death, strongly stressing the different behavior of the 
quantum world with respect to the human experience. It 
is commonly accepted that quantum behavior vanishes 
as the system size increases. Remarkable exceptions are 
quantum systems which undergo a phase transition, as 
superconductors and superfluids. A quantum superpo- 
sition of mesoscopic states has been observed in SQUID 
devices Q and seems to be a promising tool for the real- 
ization of a quantum computer. 

Quantum macroscopic states are expected to be robust 
with respect to decoherence and thus ideal candidates for 
quantum information storage. Moreover, it is well known 
that a quantum system which undergoes a phase transi- 
tion lives in one of a particular set of states, for a time 
which becomes infinitely large in the limit of large system 
size. In particular, if the ground state is twofold degen- 
erate, one can associate these states to a macroscopic 
quantum bit. The availability of macroscopic quantum 
bits is relevant for quantum information processing, as 
shown, for instance, in Ref. |3j, in the case of spin clus- 
ters. Here a new application in teleportation processes is 
also shown. 

Decoherence of a single qubit has been extensively 
studied 0, H II SIS IST3] One of the most relevant 
causes of decoherence is the coupling with a bosonic bath 
El 13 ES El 13 13 113 whose effects are relevant also 
at zero temperature. 

In the present paper we investigate the coherence of 
an array of N double quantum dots coupled through 
Coulomb interaction in order to show that such system 
is a suitable candidate as a macroscopic qubit. The first 
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step is to show that in the long time limit the array be- 
haves as a two level system with energy separation which 
vanishes for large N. The model is exactly equivalent 
to a one-dimensional antiferromagnetic Ising model in 
a transverse field. The system is characterized by two 
equivalent charge configurations and in the macroscopic 
limit the ordering at zero temperature implies a sepa- 
ration of phase space in two regions around each of the 
degenerate configurations. However, for a finite size sys- 
tem, hopping induces oscillations between the two con- 
figurations. We associate to this behavior a macroscopic 
quantum bit. The antiferromagnetic Ising model in a 
transverse field has been studied since the pioneering 
work of Bethe 0, 0, |2(| , and has received recently a 
renewed attention as a model for quantum computation 
pll |2~3 . |3|. We found convenient to introduce a sim- 
ple approximation which makes transparent how the two 
level behavior appears asymptotically. 

The study of decoherence in such a system is analo- 
gous to decoherence in a quantum register [24| . We show 
that, at least in the weak coupling and zero temperature 
limit and for a three-dimensional environment, the sys- 
tem exhibits a robustness growing with the size of the 
array. 

Decoherence with respect to phonons of a single two 
level system has been studied with various methods (see 
reviews by Leggett et al. E3 and Weiss E3)- We found 
however the resolvent method [3 H3- already intro- 
duced in the discussion of electron-phonon interaction 
problems [2]]], convenient to obtain results at zero tem- 
perature in the double dot chain. 

The paper is organized as follows. In section [H] we 
introduce the model of the double dot chain and its in- 
teraction with a phonon bath. The introduction of the 
resolvent method to discuss decoherence will be the argu- 
ment of section lTTTl In section lTvl we shall apply the same 
method to show that the double dot chain, in the limit 
w/U << 1, being w the hopping amplitude between dots 
inside each pair and U the Coulomb interaction between 
different pairs, behaves as an effective two level system 
with energy separation decreasing exponentially with N. 
In section we study decoherence in our system in the 



2 



approximation introduced above. Finally, section IVII is 
devoted to conclusions. In appendix we review, inside 
the present approximation, decoherence effects for a sin- 
gle dot pair. 



II. THE MODEL 



H s = -w £ of + \ £ (af af +1 + l) , (4) 
i i 

H B = ^W q a q (Z q , (5) 

q 

H SB = Vj q ¥" cosM K + «- q ), (6) 



In a previous work [28j we proposed an array of few 
coupled quantum dot pairs as a channel for teleporta- 
tion. We want here to show the extension of the model 
to a case of N pairs and discuss the robustness of the sys- 
tem with respect to decoherence due to interaction with 
an external phonon bath at zero temperature. We ex- 
pect that the extensive character of the interaction will 
increase decoherence while the macroscopic nature of the 
first two energy states will enhance the robustness of the 
system. It will be shown that the latter feature prevails. 
Neglecting spin effects, the double dot array is character- 
ized by the Hamiltonian 



H S = U 



N-l 2 



1=1 a=l 



ni a n l+1 



N 

i=i 



(1) 



where c\ a creates an electron on the 1 (th) dot on the a 

(th) row of the array and ni ia — c, „cj )0 . 

To extend the teleportation scheme described in |28| . 
we introduce an initial superposition of two spin config- 
urations of zero potential energy of N pairs: 



\s) = a |$)+ 



(2) 



where |$) = ||, T, I, T, •••• T) and |*) = |T, I, T, I, •••• I)- 
Let us consider an initial system of N— 1 double quantum 
dots with hopping inside any pair and without Coulomb 
interaction. The ground state of this system is repre- 
sented by the tensor product of (|f) + j\/2 for each 
pair. By an adiabatic switching of electrostatic repul- 
sion between adjacent pairs, the system is driven in its 
new ground state, which, for w/U << 1, is well approx- 
imated by (|$jv-i) + \^n-i))/V2. \$n-i) 
is the same state than |$) (I*)), defined on N — 1 
sites. The state \S) is obtained considering an extra dou- 
ble dot (as usual called Alice) in a superposition state 
a |t) + (3\l) and its interaction with the first pair of 
(|$jv-i) + |*iv-i))- If the interaction is adiabatically 
switched on again, then the system is driven in a state 
close to \S). A proper manipulation of system param- 
eters permits to transfer the information, i.e. a and /3, 
encoded previously by Alice, to the last double dot (Bob) 

The model described above is suitable to be rep- 
resented by a spin Hamiltonian through the mapping 
of = (rt;,! — rt; j2 ) and of = (cj^cj.a + h.c). This pic- 
ture is useful to study the decoherence effects induced by 
the interaction with a phonon bath. In the spin repre- 
sentation the overall Hamiltonian becomes 



H 



Hs + Hb + Hsb, 



where a mapping between the spin states ||) and |J.) 
and the charge states |1, 0) and |0, 1) has been performed 
and ni = (erf + 1) /2. We indicate with 9 the angle be- 
tween the phonon mode q and the the dot chain direc- 
tion. This notation is useful for describing a generic d- 
dimensional environment coupled with a one-dimensional 
system. The constant <? q represents the coupling of the 
dot charge with the mode q. The explicit mathematical 
expression for c/ q depends on the specific configuration of 
the system and the type of interaction. In Ref. [2^ the 
explicit form of g q in some remarkable case is given. 



III. RESOLVENT METHOD FOR WEAK 
COUPLING DECOHERENCE 

Decoherence at zero temperature is studied using the 
resolvent method. At the initial time t = system and 
bath are decoupled: |H (* = 0)) = \S) ® |0) where |0) is 
the vacuum phonon state. 

The time evolution of the state |3 (t)) = 
exp (— iHt) |H (t — 0)) is studied in terms of the 
complex Laplace transform defined as 



|3 (u>)) — i lim 



Auit—ot I" 



3 (*)) dt. 



The resolvent method allows to write 

1 



\S(uj)) = 



to — H 



Using the identity 
1 



H 



H 



[H(* = 0)> 



H 



H 



(7) 



(8) 



(9) 



and performing a projection on the vacuum phonon state, 
we define a new system state |$s (u)) = (0|H(cj)) that 
obeys to the evolution equation 



|$ s («)) = 
1 



1 



uj - H s 



-H. 



SB 



uj- H 



\®s (* = 0)) + 
|S(* = 0)). 



(10) 



Here the bath ground state energy is set to zero and 
H = H S + H B and Hi = H SB . 

In the weak coupling limit only corrections to the 
imaginary part of |$g (u>)) will be taken into account. 
We first perform an iteration inside Eq. 1101 replacing 
(uj — H)~ with the right hand side of Eq. ^\ and then 



(3) introduce a complete set of intermediate phonon states: 



3 



l*sM) = - J 7rl*s(t = o)) + (o|^— - 

to — il5 to — -H5 



H. 



1 



si; 



|S(* = 0)) 



i 



to - Hs 



-H. 



1 



SB 



w — Hs — Hb 
I 



uj — Hs — Hb 

H SB \k) (k\^—\S(t = Q)) 
uj — H 



(11) 



In a perturbative approach, terms involving powers of g q 
are small and can be neglected, unless self-energy con- 
tributes appear. In the latter case, a not negligible imag- 
inary part can arise performing the sum over q in the 
continuous limit. 

The contribution involving self-energy in the sum cor- 
responds to k = 0, since (0| (to — H)~ |3 (t = 0)) is ex- 
actly |$s (to)). Then, all other linear and quadratic con- 
tributions in Hsb will be neglected. Hence, Eq. 1111 be- 
comes 



1 



1 



uj-Hc 



G(H S ) \f> s M) = 



1 



uj - H f 



\*s (t = 0)) , 
(12) 



where 



G(H S ) = (Q\H SB - 



1 



Hs-Hb 



H sb |0) 



(13) 



is the self-energy operator acting on the system sub- 
space. The right term of EqEI describes the evolu- 
tion of the macroscopic state isolated from phonons. As 
we will show in section llVl in the limit of w/U <C 1, 
the macroscopic dot chain behaves as a two level sys- 
tem oscillating between the Hs asymptotic eigenstates 
|±) = 2- 1 / 2 (|$) ± |*)) with energies E±. So, Eqdbe- 
comes 



1 



uj- H s 



G(H S ) |$ s M> 



uj-E 4 



|+) (+|4 5 (t = 0)) 



uj - E- 



-)<-|*s (* = 0)) 



(14) 



Noting that the operator G(Hs) maps the subspace 
spanned by |±) into itself, it is possible to reduce EaH4l 
in terms of two coupled equations: 

{uj-E+- C++) <+|$ s (uj)) - G+- (-|$ s {uj)) = 

(+|<M* = 0)>, (15) 

(uj - E_ - G ) (-|$ s (uj)) - G-+ (+|$ s (w)> = 

(-|* s (t = 0)>, (16) 

where G ±± = (±|G|±). 

To the leading order in the system-bath coupling, we 
obtain 



(+|$sH) = 

(-|$s (uj)) = 



1 



uj - E+ - G^ 
1 

u> — E- — G~ 



(+\$g(t = 0)), (17) 
(~|<i>5 (t = 0)).(18) 



The solution in the time domain is obtained assuming 
first the correction introduced by the matrix elements 
of G as negligible, and then calculating the latter in in 
uj = E+ or uj = E- . 

For instance, the integral 



-duj 



E + - G++ 



is calculted assuming first G ++ = 0, obtaining for the 
pole uj = E + , and then substituting this value inside 
G ++ , which depends on uj. After, the principal value of 
G ++ will be ignored, and only the imaginary part will 
matter. We compared, in appendix, the results of our 
approximation, with those known for single quantum dot 
pairs. 



IV. DOUBLE DOT ARRAY EVOLUTION 

As first step we calculate the evolution of the system 
when it is decoupled from the bath. We find convenient 
to explicitly solve the evolution from an initial state cor- 
responding respectively to |$) or |^) introduced in Eq. 

m 

We distinguish in the system Hamiltonian the hop- 
ping term Hi = —w^2icrf from the potential energy 
H = U/2 (vf^i+i + l)] • This is the antiferromag- 
netic version of the well known one-dimensional Ising 
Model in a transverse field [3l]]. It is worth noting that 
the absence of periodic boundary conditions implies a 
relaxation mechanism of an initially ordered state where 
a single domain wall propagates between the two end 
points of the array. This feature makes a difference in 
the excitation spectrum which is relevant for an array of 



4 



finite size. 

Applying Hi on |$ (t = 0)) the system is driven in a 
new configuration labeled as |$i (t = 0)). The action of 
Hi generates a sum of states each of them differentiates 
from |$ (t — 0)) due to one spin flip in a different place 
along the array. Here it is important to note that flips 
on the first and the last qubit put the system in a state 
with Coulomb energy U , while all intermediate transi- 
tions lead to a state with a 2U electrostatic energy. In 
the limit of U large with respect to w, we shall neglect all 
configurations involving intermediate states with energy 



greater than U. 

In each step of a repeated application of Hi it is pos- 
sible to go towards new configurations or to come back. 
Then, for n > 0, we write 

Hi |$„ (t = 0)) = -w [|$„_x (t = 0)) + |$„+i (t = 0))] . 

(19) 

After N steps the system reaches |\E') and after 2N 
steps it comes back to the initial configuration. Defin- 
ing |$jv) = and |$ ) = |^2jv) = |$), and taking into 
account the time evolution we obtain 



(w - U) |*„ M) = |*„ (t = 0))-w [|$ n _i (a;)) + |$„+i (w))] - C/ (£„,o + tfn.jv) |*„ (w)) 

I 



(20) 



The system is solved by means of the discrete Fourier 
transform defined as 



2JV-1 



The asymptotic behavior is determined by values of ui 
close to zero. Then a(uj) << I and the denominator of 
B n (w) reads as geometric series: 



V n— 

2JV-1 

!*«(«)) = ^ E 



ink 



fc=0 



I [/ 



2JV-1 



e m » 9 ^ a' (w) cos J -j^g, 



-ink 



(21) 



g=0 



(27) 



or 



As a consequence of periodicity conditions, k = be- 
ing n = 0, 1, 2, 2iV - 1. 
From Eq. [3U| follows 



B.M4A™EtE(:)(f' 



2Nuj-U 



[oj - U + 2w cos fc] $ fe (o;)\ = $ fc (t = 0) 
J7 



9 =0 /=0 m=0 

x exp |i— (2 — 2m — n)q . (28) 



/2iV 



(j$o(w)) +e iArfc |$Ar(w))) . (22) The sum over g gives 



It's now possible to extract two equations connecting 
|$o) to |$jv): 



|$ H): 

where 



[1 + B (to)} \A (w)) - Bjv (w) |Av M) 



(23) 



" M ~ 2N lo - U ^ ^ [ m J 2 l 1 - e ** ('-»»-») 
;=o m=0 v 7 

The condition for a nonvanishing _B„ (u>) is 
(2 — 2m — n) = 2NK, where K is any integer be- 
tween — oo and +oo: 



[1 + B M] \A N (w)) - Sjy (w) | A) (w)) 



(1 + Bo H) 2 -^ ( w ) 



I A, M> = 



E 



$ fe (t = 0) 



/ 2A r ^ uj -U + 2wcosk 

h—Q 



(24) 



(25) 



[/ « a l (u,) 

Bn M Z^U ^ ^ m\(l-m)\ 2' 5 ('-*»-"> 
;=o m=0 v y 



or, using the Kronecker Delta function 



,2NK- 

(30) 



S„ (w) 



I U 

2Ntv- U 



2N-1 



y — - 

^ l-a 



9=0 



(w) cos J^g 



(26) 



B n (w) — jj ( l+n+2NK \ \ ( l-n-2NK \\ o' ' 

1=0 K=-oo V 2 7- V 2 r 

(31) 

Since the coefficients of a Newton's binomial formula have 
to be real and positive, in the limit a (u>) << I we obtain 



with a (uj) = 2w I (U — oj) and noting that Bn = B- 



B (to) 



U 



oj-U 



(I + M) . 



(32) 



5 



where 



2JV-1 



M 



1 ^ 1 
1 ~ 2N 1^ 1 - ^ 



q=0 



cosq 



(33) 



contains powers of w/U and has to be calculated at the 
desired order in q, and 



N 



(34) 



Here we note that the last contribution cannot be ignored 
because it gives rise to the energy separation between 
|$ ) and |$jv). 
Furthermore we obtain 



and 



\A ( W ))~i|* (« = 0)) 



l^iv (t = 0)> 



(35) 



(36) 



As a result, after an inverse Laplace transform, we get, 
apart from corrections containing powers of w/U, 



l*o (*)> 



AMUt 



|$o (0))cosAi-i |$jv (0)) sin At] 



and 



(37) 



|$jv (f)) = e lMUt [|$jv (0)) cos At - * |$ (0)) sin At] , 

(38) 

having introduced the energy gap 



A = 2w{2w/U) 



N-l 



(39) 



We eventually obtain the long time behavior of a two 
level system with energy separation exponentially van- 
ishing in the large N limit. Actually, in Ref. [ijj (see 



equation (3.32c)) the eigenvalue of Eq. 1391 was derived. 
On the basis of this result the phenomenon of asymptotic 
degeneracy was established and shown to be directly re- 
lated to the appearance of the ordered phase in the large 
TV limit. 



V. DOUBLE DOT ARRAY DECOHERENCE IN 
THE LONG TIME LIMIT 

According to the previous analysis we can limit ourself 
to consider only the first two states |±) of the array. 
The decoherence rate will be however modified by the 
extensive interaction with the bath. 

We have to calculate the matrix elements of G {Hg) 
in the subspace of |+) and |— ) taking into account the 
particular system-bath interaction Hsb defined in Eq. 
Here 



G{Hs) = J2 ei9C ° Se{l ~ l,) \^\ 2 nv 



1 



LO — Hs — to, 



-ni, 



(40) 

where the sum over 1,1' runs over the array sites where 
electrons are present. 

We choose the basis elements |+) and |— ) defined re- 
spectively as the sum and the difference of |<£>) and 

We introduce the form factor A q cos g defined through 



J2n l e i i cosei \$) = A qcoB g\$) 



(41) 



or 



J2n l e i « cos6l \y} = e i 0A qcos g\y) (42) 



Explicitly, A qcos9 = (l - e l2qcosf)N ) / (l - e i2qcosf) ). 
The matrix elements of the self-energy operator are 



J 



q 

G = £|<? q | 2 |A 



2 i a |2 

cos 



_,2 q cos 9 



2 q cos ( 



LO — E + — UJq LO — E— — LOq 



q cos f 



COS 



2 q cos t 



2 q cos 6 



LO — E- — LO a LO — E+ — LO c 



q cos a . q cos f 



G + = (G+ )* = l5 q | 2 |A gC ose| 2 cos — - — sin 



LO — E + — LOq LO — E_ — LOq 



(43) 
(44) 
(45) 



We find convenient the introduction of the following generalized densities of states 

+- r\ I -+^* i \ -Sr-\ |2 1 * i 2 qcos9 . qcosO 
P ( e ) = (P ) ( e ) = -*2^l5ql l A ?cose| cos — - — sin — - — d(e-w q ) 



Pi (e) = Ifql 2 l A gcose| 2 cos : 



gcosf 



(46) 
(47) 
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P2 (e) 



^2 l5q| 2 l A 9co SI 
q 



2 . 2 gCOsg 

I sin — — <5(e-w q ) 



from which follows 



G+- 
G— 



= / 



Pi (e) 



P2 (e) 



= / rfe 



/ / dep+- (e) 
Pi (e) 



e uj — E- — e 
1 



to — E. 



- 



e uj 



P2(e) 



(48) 

(49) 
(50) 
(51) 



The real part of G gives a negligible contribution to the 
pole location if compared with E_ and E + . Thus, assum- 
ing a density of state different from zero only for positive 
e, as m Ref. 0, the only non vanishing contribution is 
7= -ImG— : 

7 = -7rp2 (A) , (52) 

where A = £L — E+ is the energy gap of the two level 
system and is positive (being |+) the ground state). 
Then the solution for (+|$s (*)) and (-|$s (t)) is 

(+|<M*)> = e iE +*<+|$ s (t = 0)) ! (53) 
(-|* s H) = e lB - t -^ t (-|$ s (< = 0)). (54) 

As expected, the ground state is not affected by de- 
coherence, while the excited state relaxes. Damping is 
proportional to the density of states calculated at the 
energy gap. The density of states is however quite dif- 
ferent from that of a single dot pair. Two competitive 
effects appear. The first one is represented by the pres- 
ence of the form factor A gcos e inside p2, which, in the 
large N limit, increases the dephasing rate by a factor 
proportional to TV 2 . The second, predominant, effect to 
be considered is the expontential reduction with N of the 
energy separation. 

For instance, in the simple case of |g q | 2 = 1/N and 
w q = cq (longitudinal phonons) 

7(A)tx d cos 6d d q — ^— sin 2 ^— — <5 (A - c q ) , 

7 sin g cos 2 

(55) 

where d is the dimension of the bath and c is the speed 
of sound. If we compare this quantity with the system 
oscillation frequency we obtain 

^ cx AT 2 A*/ 2 " 1 . (56) 

This result indicates that, for a phonon bath in three 
dimensions, the macroscopic limit involves a growth of 
the robustness with respect to decoherence. 



VI. CONCLUSIONS 

The existence of a macroscopic degenerate ground 
state is a general feature of a system exhibiting a phase 
transition. It would be of interest to exploit such degen- 
erate states as elements of a macroscopic qubit. Here 
we considered an array of N interacting dot pairs. We 
showed that, in the long time limit and neglecting cor- 
rection of order w/U, the system oscillates between the 
two configurations characterized by zero electrostatic en- 
ergy. Such a system has a robustness which increases 
with the size as shown calculating the decoherence due 
to a phonon bath at zero temperature. Decoherence cal- 
culation has been performed in the framework of an ap- 
plication of the resolvent method. It is however impor- 
tant to note that a general conclusion about robustness 
of macroscopic qubits should consider nonzero tempera- 
ture effects and other different mechanisms for daphasing 
in semiconductors (such as cotunneling and background 
charge fluctuations). 

As said, coherent manipulation of \S) is useful in order 
to realize a support for quantum information transfer al- 
lowing a solid state teleportation. The basic information 
processing steps, i.e. initialization of the system, local 
gates and readout are respectively obtained by adiabatic 
variation of system parameters, oscillations between |$) 
and 1^)) and local charge measurement. Moreover, it is 
worth to note that, because of the analogy with spin clus- 
ters behavior, most of considerations on quantum com- 
puting developed in Ref. Q should be valid also for our 
system, with the advantage that instead of measuring 
spin states, we need to detect local charges. Due to the 
asymptotic behavior of the effective hopping amplitude, 
such a qubit requires gate times exponentially increasing 
with the size of the system and the optimal chain length 
will be determined by the specific application. 



APPENDIX A: DECOHERENCE RATE IN A 
DOUBLE QUANTUM DOT 

We introduce a double quantum dot in contact with a 
bosonic bath with Hamiltonian 



7 



H — Hs + Hb+Hsb-, 
H s = |cr 2 + T(J X , 



H 



h sb = \az^9q{a\ 



(A-l) 



discussed in Ref. [12|. Here a one-dimensional bath is 
considered for simplicity. Labeling with \L) and \R) the 
eigenstates of a z with respective eigenvalues +1 and— 1, 
the eigenstates of Hs are 



\±) = —{±2T\L) + (A T e)\R}}, (A-2) 



where A = Ve 2 + 4T 2 and A± = ^(A^e) 2 + 4T 2 while 
the respective eigenvalues are e± = ±|A. 
By inversion we obtain 



2V, 



A 



4TA 



A_ 



4TA 



(A-3) 



Iff} - N+ 
lR) 2A 



A_ 
2A 



(A-4) 



Eq. E| has now to be solved using 



G(Hs) = lY,\9 q \ 2 °z 



(A-5) 



We need to calculate (+| G (H s ) |+) and (-| G (H s ) |-). 
Actually, obtaining G ++ will be enou gh, due the intrin- 
sic robustness of the ground state |— ) |3(| which implies 
that G~~ has to be zero (this feature is easily checked in 
the present formalism). To do it first we write |+) in the 
\L,R) basis, then apply er z , come back in the |±) basis 
in order to apply (to — uj q — Hs) -1 , rewrite the new state 
through |L, R) to apply the second a z operator, and fi- 
nally re-express the result in terms of |+) and |— ). The 
result is 



= i Em 2 *_a(U 



bJ — UJ„ + 



#( A A 



(A-6) 



The sum over q is performed as an integral through the 
introduction of the density of states p which is assumed 
to be different from zero only for positive values of its 
argument . The second term inside the square bracket 
gives the contribution to the imaginary part, which is 
7 = - [nT 2 p (A)] /A 2 . The evolution is thus 



The density matrix in the basis |±) is then 



PW - \ p -+ (o) e -^e jAt 1 - p++ (0) e- 2 ^ > ' [A ~ y) 



(+|*s (t)} = (+|$ s (i = 0)) e-^e"'* 



■p(A)t 



1*5 (*)> 



-\*s (t = 0))e l 



(A-7) 



(A- 



with the same dephasing rate obtained in |T^. in the 
regime of zero temperature, using markovian assump- 
tions. 
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